UCLA ENG-8746 
December 1987 


UTILIZATION OF FEP ENERGETICS 
FINAL REPORT 
December 31, 1987 


by 


T. H. K. Frederidng, Principal Investigator 
P. Abbassi 
F. Afifi 

P.K. Khandhar 
D. Y. Ono 
W. E. W. Chen 


/)/?) i 

//J 



PREPARED FOR THE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
AMES RESEARCH CENTER 
MOFFETT HELD, CA 94035 

Grant NAG 2-464 


Department of Chemical Engineering 
School of Engineering and Applied Science 
University of California, Los Angeles, CA 90024 


(NASA-CR-182373) UTILIZATION OF FEP 
ENERGETICS Final Report {California Uni*.) 

P CSCL 1 3 K 


G3/37 


N88-14379 


Unclas 

0118109 


UCLA ENG-8746 
December 1987 

UTILIZATION OF FEP ENERGETICS 
FINAL REPORT 
December 31, 1987 


by 


T. H. K. Frederking, Principal Investigator 
P. Abbassi 
F. Afifi 

P.K. Khandhar 
D. Y. Ono 
W. E. W. Chen 


WITH APPENDIX SECTIONS BY S. Caspi, S.W.K. Yuan et al. 


PREPARED FOR THE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
AMES RESEARCH CENTER 
MOFFETT FIELD, CA 94035 

Grant NAG 2-464 


Department of Chemical Engineering 
School of Engineering and Applied Science 
University of California, Los Angeles, CA 90024 



TABLE OF CONTENTS 


1. ABSTRACT 

2. INTRODUCTION 

3. FEP SYSTEM 

4. CYCLE STUDIES 

5. POROUS MEDIA TRANSPORT CHARACTERISTICS 

6. SUMMARY INCLUDING COUNTERFLOW COOLING OF COMPONENTS 

7. REFERENCES 


PAGE 

1 

2 

6 

14 

25 

26 
38 


APP. A: STABILITY AND REFRIGERATION OF MAGNET CRYOSYSTEMS 
NEAR 1.8 K USING THE THERMOMECHANICAL EFFECT 39 

APP. B: VORTEX SHEDDING REGIME IN POROUS MEDIA 

APP. C: EXPERIMENTAL DATA (POROUS PLUG DATA) 

APP. D: PERMEABILITY DETERMINATION 

APP. E: THERMOMETER CALIBRATION AT VERY SMALL TEMPERATURE 
DIFFERENCES 

APP. F: ABSTRACTS OF RECENTS CONFERENCES 

APP. G: REFRIGERATION AND HEAT PUMP SYSTEMS BASED ON HE II 
VORTEX CONTROL 

APP. H: SPACE CRYOGENICS COMPONENTS BASED ON THE 
THERMOMECHANICAL (TM) EFFECT 

APP. I: DATA OF THERMOMECHANICAL HEAT TRANSFER EXPERIMENT 75 


43 

45 

52 

53 
55 

57 

68 



1 . 


ABSTRACT 


The R&D work on Fountain Effect Pump systems (FEP systems) has been of interest 
in the competition between mechanical pumps for He II and FEP units. The latter do not 
have moving parts. In the course of the work, the energetics has been addressed using one 
part of a simple "4-changes of state" cycle. One option is the FEP ideal change of state at 
constant chemical potential (p). The other option is the two-state sequence p-P with a dp = 
0 state change followed by an isobar. Questions of pump behavior, of flow rate response to 
temperature difference at the "hot" end, and related questions of thermodynamic cycle com- 
pletion and heat transfer have been addressed. Porous media data obtained elucidate 
differences between vapor-liquid phase separation (VLPS) and Zero Net Mass Transfer 
(ZNMF). 
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INTRODUCTION 


The project has been started up with the goal of quantification of the energetics of 
fountain effect pump systems (FEP systems) which utilize the thermomechanical effect in 
liquid He II in the superfluid range. Mostly 1.2 to 2.1 K is a rough characterization of tem- 
peratures of interest in various large scale systems. In the course of the work it became clear 
that the topics ought to include both the equilibrium state thermodynamics and the transport 
phenomena of the process components. The emphasis has been on micro-gravity noting that 
the first explicit successful use of thermomechanical forces (TM forces) has been in space 
vessels in vapor liquid phase separators (VLPS). The IRAS system has become famous in 
this respect The ongoing efforts focusing on additional controlled fluid management reach 
out farther. An example is the transfer of liquid in space using the FEP type systems without 
moving parts. Prior to a detailed inspection it appears to be useful to focus on historical 
developments. 

1988 may be considered the half century jubilee year for superfluidity research and 
development. Though indications, such as the lambda transition in liquid Helium-4, have 
been known since Dana and Onnes work, the systematic study of superfluidity did not start 
before about 1938. In that year the thermomechanical effect had been reported by Allen and 
Jones 1 . Nearly simultaneously, superfluidity was found, i.e. the ability of He EL to flow 
through very narrow channels. These findings had been the subject of sort of a race as both 
Kapitza 2 and Allen and Misener 3 got data in the area of inquiry. Concomitantly the two-fluid 
model came along with a surprising speed in view of entirely new phenomena, the 
"supereffects" of He n. In particular, the London version became accessible to Western 
researchers while the Landau version remained in the Eastern domain at first. Thus, major 
theoretical developments gave immediate support to data interpretation efforts. 


2 



It did not take long time, after the advent of the Collins liquefier, to launch demonstra- 
tions of the fountain effect, and of the static version of the thermomechanical effect, in major 
universities having "Collins'’ systems. Thus, many additional phenomena of the underlying 
physics got significant attention. An example is the resolution of first sound in He II, the 
usual pressure density disturbance propagation, and of second sound, the temperature entropy 
mode. Other sounds modes followed e.g. third sound, fourth sound and fifth sound. 

A comprehensive program of research had been started by Hammel and Keller 4 in the 
area of He II film flow and in thermodynamics and transport phenomena. In particular very 
fine slits were used to see the fountain effect thermostatics. According to London’s equation, 
the thermomechanical pressure increment is given by 

dP = p S dT (2.1) 

S = entropy per unit mass, p density. 

Therefore, the difference in pressure due to a temperature difference A T is 

AP = Jp S dT (2.2) 

In the context of this equation, Hammel and Keller obtained limits imposed on the maximum 
value of AP attainable. Also critical velocities appeared to be an obstacle toward immediate 
full utilization of the pressure difference predicted by equation (2). Indeed kinetic energy 
effects had to be taken into account. An example used in the two-fluid frame of equations is 

V(p + v&) = 0 (2.3) 

The kinetic energy, per unit mass, of flowing superfluid is v s ^ 2 ; p is the chemical potential per 
unit mass, v s superfluid velocity. The entire set of data has not been easily accessible to sim- 
ple theory as there were numerous critical velocity values for various modes. In this context, 
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the entropy was checked in detail because of concerns of entropy contributions not assigned to 
normal fluid alone. Nevertheless, a special mode called zero net mass flow was among the 
early systems studied. It turned out that He II had an enormous axial heat transport capability 
despite operation in the "supercritical" regime, i.e. beyond the critical velocity. Thus, the 
thermomechanical forces (TM forces) appeared to be at work, and interesting arrangements 
were created to take advantage of this type of "heat piping" properties of He n. Thus, in a 
sense, the "first use" of VLPS, as an application of He II TM forces, has been superseded in a 
subtle manner by the laboratory use of zero net mass flow conditions. 

Coming back to devices, it is noted that F. London 5 already looked at thermodynamic 
state changes of a special system. One may consider this a "thought experiment". One may 
also look at this as a first proposal for a system with the potential of practical realizability. In 
1938 things were very far away from such a status. In fact occasionally the opinion had been 
voiced that Equation (2) does not provide a significant incentive for practical use. The argu- 
ment goes as follows: The critical velocities are mostly very low; the pressure difference 
attainable at the critical velocity is not large; further, the energy input in the form of heat does 
not provide for a "very efficient" system from the second law of thermodynamics point of 
view. These arguments contain some truth however they ignore several extremely interesting 
features for the space fluid management. 

First of all, the low AP and low critical velocity may be overcome by suitable very 
fine porous plugs and in particular by multi-stage pump systems. Such a system did not come 
up before the Eighties. Severijns 6 demonstrated the use of a 4-stage pump which had very 
good agreement with predictions based on Equations (1) and (2). The second major point is 
the option of using "waste heat" as input, for instance parasitic heat inputs available through 
most cryogenic insulation systems e.g. in supports, current leads and others. The second law 
of thermodynamics is valid. However it provides more of a comparative measure for various 
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schemes, Considered competitors, than a serious limitation on the "waste heat" use. Last not 
least, the question of realiability, by elimination of as many moving parts as possible, is 
another factor. The FEP unit permits pressure built up without moving parts. Thus it has an 
essential advantage over other options. 

It is possible in principle to push device frontiers solely by many empirical develop- 
ments emphasizing components of interest However, theoretical guidance is valuable for a 
reduction of the total numbers of data sets required. 

An early example is the laminar heat flow mode at Zero Net Mass Flow (ZNMF). 

Tisza 7 argued in 1947 for a heat flow rate q proportional to the mean speed of normal fluid, 

proportional to the pressure gradient, Eq. (2.1), and to pST and (ps/pn); ps is the superfluid 
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density, pn normal fluid density. The work of Keller and Hammel and of other researchers 
showed that (ps/pn) does not enter. It influences q however in the "supercritical" regime. 

A more recent example of this nature has been VLPS work. Yuan 8 clarified important 
controversial points. The rate - controlling mode is flow of heat in the form of normal fluid 
flow in porous media separators. The recent book of Wilks and Betts 9 has documented pro- 
gress in the field. Also it lists simple two-fluid equations referred to subsequently. The 
present report contains various sections as given in the table of contents. Appendix sections 
listed in the Table of Contents describe further experimental details and additional studies 
conducted in the present context. 
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3. FEP SYSTEM: 


PUMP. PLUG STUDIES WITH SERIES CONNECTION OF PACKED A1 2 0 3 PLUG AND 
Cu POWDER HEATER 

There are two items: 

a. Flow rate dependence on driving temperature difference; 

b. Temperature profiles of the combined packed powder-heater system. (FEP runs) 

The flow rate dependence work dates back to the days of near-isothermal flow investigations 
of Vote et al 10 . The thermomechanical rate equation has been adopted as a modified power 
law equation which contains recent findings concerning an "effective theimomechanical gra- 
dient". This latter version is a result of studies by Dr. Sidney Yuan et al 11 . For easy refer- 
ence only the simplest cases are included here giving an example of the "1/4 power law" in 
each case. 

The original Vote et al 10 , work did not refer to the permeability of a particular flow 
system geometry. Instead the hydraulic diameter (D h ) of slits and ducts with circular cross 
section has been used. The equation is given as superfluid mean speed v s , or as mass flow 
rate m. 

Superfluid mean speed 

v s = const [ I AP I /D h ] 1/4 p -3/4 ri n 1/2 (3.1) 

Mass flow rate: 

m = A c p s Const [ I AP I /D h ] 1/4 ri n 1/2 /p 3/4 (3.2) 

(A c cross section; p density, VP Pressure gradient for flow induced by external pressure 
difference or by gravity forces with difference AP = gp AZ; Z position). The equations 
appear to be "guide posts" for other cases if a switching to the thermomechanical P-gradient, 
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compare (Eq. 3.1) would suffice. Another convenient option appears to be the use of L c = 
K 1/2 ; (k Darcy permeability). L c replaces the hydraulic diameter. It is useful in particular for 
porous media. 

An objection to the Vote et al approach is the observation that for several flow cases it 
is the "effective thermomechanical gradient" (VP T ) (ps/pn) which governs VT - induced 
flow. 

a. Flow rate dependence on driving temperature difference 

Any introduction of thermal quantities will in general cause departures from the near- 
isothermal resistance evolution described as power law by the Vote et al equations. The latter 
have, at any point of steady operation, a well-defined flow resistance. If the resistance is 
above the normal resistance or equal to it, the nature of the flow is Newtonian. There is 
either laminar flow or turbulent flow. The latter usually is roughness dependent. If the flow 
resistance is below the classical normal value, the system is operating on a point of the resis- 
tance transition from the superfluid to the normal flow states. 

As soon as an externally applied temperature difference is switched on, the system is 
capable of responding with a flow rate which is not necessarily characterized by the simple 
arguments presented above. There are not enough data available at present to clarify this 
aspect of superfluid mass flow. However, the equations of Vote et al. have been modified by 
Yuan in order to characterize various FEP plugs. The equation incorporates the effective 
driving thermomechanical gradient (Ps/p n ) p S VT. The superfluid speed and the mass flow 
rate are given as power law quantities for the "1/4" law as analogs of the preceding equations. 
Other power laws arc referred to subsequently. 

Superfluid velocity (mean value superficial): 

v s = K 4 (I A P t I/L c ) ,/4 (p s / p n ) I/4 p" 3/4 Ti n ,/2 (3.3) 
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or 


(3.4) 


* = V s p, A c - A c p s K„ ( I A P T I /LJ 1 ' 4 (p, / Ph ) 1 ' 4 p-“ r | n 1/2 

The quantity K 4 is a constant for each pump plug as long as the power law is valid, i.e. for 
moderate vortex shedding rates. L c is the square root of the Darcy permeability. Insertion of 
the effective gradient (eq. 3.4) permits a form of the equation quite analogous to the preceding 
Vote et al equation (3.1). 

Mass flow rate 


m = p v A c = p s v s A c =* K 4 [IV P etf l]" 4 TV 1 ' 2 p- 3 ' 4 p s A c (3.5) 

The preceding equation is to be used noting that pv = p s v s , and that A c is the total plug cross 
section. Pump plug data for K 4 appear to be much less dependent on the permeability than in 
VLPS. The pump plug results for K 4 given by Yuan et al 11 are on the order of unity for per- 
meabilities below 10 -10 cm 2 . This implies that they are very much larger than corresponding 
VLPS plug constants for the same permeability. 

Other power law exponents have been quoted in the literature. The most frequent 
number is d log m/d log AT = 1/3 ("derived" from the Gorter Mellink mutual friction force 
for zero net mass flow). For the near-isothermal transport of the Vote et al. equations, things 
are quite readily distinguished. The exponent d log RESISTANCE/ d log AP varies from 
unity via 3/4 to 2/3 monotonically with an increase in flow resistance. When n is converted 
to the exponent for mass flow versus pressure difference, we obtain m = d log v/d log AP = 1 
- n. For instance, for n = 2/3, m will be 1/3. In contrast, the classical Blasius value is 
characterized by n = 4/7 = 0.5714. It is noted that the classical friction factor (f) is character- 
ized by d log f / d log v = - 1/4. Further, d log AP / d log v is equal to 2 - 1/4 = 1.75 = 7/4. 
The value of m is 5/7 = 0.714. In other words, the classical range of power law exponents is 
very incompatible with the quantum liquid resistance transition under consideration. 
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For the thermomechanical pump plug, again the exponent m = 1/3 has been quoted. 
There are not enough data available for consistency checks. However, if things are similar to 
the resistance transition described by Vote et al., the exponent 1/3 ought to be observed for 
high vortex shedding rates and when the size dependence disappears as one goes from one 
plug to the next one in a certain pore size range. Despite the lack of data, some insight is 
provided by the subsequent discussion of a comparison of VLPS with ZNMF (section 5). 
This section 5 has the purpose of extending the plug data base into unknown regimes. 

Further, for liquid transfer schemes at low AP, the Appendix E contains the abstract 
related to a paper presented at the Cryog./Eng. Conf. St. Charles IL, June 1987. 

b. Temperature Profiles of the Combined Packed Powder - Heater System 

This subsection reports temperature profiles of the apparatus sketched in Fig. 3.1. The 
A1 2 0 3 portion is the FEP. It accomplishes the pressure increase associated with the applica- 
tion of a temperature difference AT. Liquid He II enters from the top flowing down into the 
loop and leaving at the outlet as a fountain after passing through the J-shaped stainless steel 
tube downstream. The packed copper section is the heater section. 

The FEP tube is inside a stainless steel tube of 0.238 O.D. and 0.254 cm wall thick- 
ness. The powder is "Aluminum oxide C" of Degussa with 99.6 % AI2O3 and a BET surface 
area of 100 ± 15 m 2 /g. The average primary particle size is 20 nm. The powder has been 
dried overnight in a vacuum oven at 125°C before compression. A stainless steel sleeve and 
hardened wire rams were used to compress the powder in a step-by-step procedure involving 
15 compression steps. The void fraction has been about 50 % based on the density of 3.3 
g/cm 3 for solid A1 2 0 3 . 
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The copper powder system has been pressed into an OFHC copper tube (O.D. 0.635 
cm, I.D. 0.335 cm, length 2.2. cm. Heater windings were made of teflon-insulated constantan 
wire wound bifilarly around the tube leaving room for a carbon resistance thermometer in the 
middle of the windings. This thermometer was mounted in a slit of 0.1 cm depth. The pur- 
pose of the Cu powder was the attainment of a large surface area permitting good contact 
with fluid. The powder purity is 99.2 %. The grain size distribution is listed below for 
different mesh sizes. 


Mesh (p.m) 

Per cent 

-325 (44) 

61.6 

+325 (44) 

28.5 

+200 (74) 

9.7 

+150 (105) 

0.2 


In a step-by-step procedure 10 layers have been applied by a press in order to achieve a rela- 
tively uniform powder packing with a void fraction of about 30%. The system may be 
characterized by geometry-and heat transfer-related parameters, such as the fin length based 
on the Kapitza heat transfer coefficient h for copper. The characteristic fin length of a simple 
pin fin is given by the product of the square root of cross section A c to circumference C c and 
the square root of the length k/h. The thermal conductivity is k. The fin length is written as 
Lp = (K A c /h C c ) 1/2 . This length has been determined as 0.06 cm at 1.8 K. The heat 
transfer coefficient rises with T 3 . Therefore near the lambda point, the fin length is smaller, 
and below 1.8 K the length will become larger. Thus, there is the option of leaving a central 
empty core for enhanced flow, or a different coarse powder in the center for 1.8 K operation 
and for higher temperatures. 


11 


Temperature records are given in the next figure . They all show a flat portion and a 
subsequent decrease toward the bath temperature. The latter indicates "leakage" of thermal 
energy not utilized entirely for the thermomechanics of the FEP. Thus, termination plugs 
and/or different conditions of operation may be chosen. 

★ 

A data sample is listed below, and further details are given in the M.S. thesis of Wil- 
liam Chen 12 , 1987. 


* App. I . 
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4. CYCLE STUDIES 


The present section addresses cycles and parts of cyclic operation. The individual 
components in general ought to match the cycle option chosen for the entire system. In the 
case of the FEP, there may be a need for an aftercooler, or the liquid may be at a higher tem- 
perature without detrimental effects. For instance, in the delivery process involving a "hot" 
object to be cooled by He I, with associated vapor formation, there may be no need to cool 
fluid back to the inlet temperature of the FEP. 

In several cases, the FEP may be used as "power unit" to energize other subsystems. 
An example is the "vortex refrigerator" (Fig 4.1) using the thermomechanics and the 
mechano-calorics of the He II "supereffects". The interesting point is a lack of performance 
figures in the literature for such a refrigeration cycle. This lack indicates a deficiency in the 
matching of components for specific tasks. Another point is the recognition of the role of the 
flow mode, e.g. resistance transition, laminar flow, roughness-dependent turbulent flow and 
others. 

Kittel 13 quotes the "inconsistency", of the supercritical FEP , "with the usual two-fluid model". 
Wilks and Betts in their recent updating of the "Introduction to Liquid Helium 9 " do not go 
into particular details of the mechano-caloric effect. However Wilks 14 mentions in his mono- 
graph on "Liquid and Solid Helium" of 1967 that "even in the narrowest slits, the flow of 
superfluid will be accompanied by a flow of normal fluid". This makes it hard to accept 
"ideal" cases without knowing the real system behavior. 

Another feature relates to vortex refrigeration based on He II thermomechanics and 
mechano-calorics. Tilley and Tilley 15 consider this type of refrigerator in their first edition. 
However the second edition of their book omits the refrigerator. One may see this as an indi- 
cation of a lack of interest or of knowledge or both. Most likely, very few cryocooler special- 
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ists have taken notice of the He II vortex refrigeration option. 

Once thermodynamic changes of state are introduced as ideal changes, there is always 
a question about non-ideal conditions and tayloring of optimized components for a special 
task. As usual there appear to be different options, and each task has certain constraints 
implying a lack of generality when applied to another problem. Therefore, the following 
simplified statements are to be considered with some caution as experimental data are quite 
scarce. With these precautions in mind, cycle results are given in simplified form for the 
ideal (i-P-)i-P cycle and for the corresponding cycle with counterflow heat exchanger prior to 
the refrigeration stage. The system is shown schematically in Fig. 4.1. 

The cycle designation is based on constant parameter references. The initial state 
change dp. = 0, or \i = const, designates constant chemical potential. The pressure is raised in 
an ideal manner, as described by London’s equation. In the second state change, isobaric 
cooling is achieved. The third state change implies constant chemical potential with a depres- 
surization. This state change leads to cooling. Ideally, the decrease in temperature is given 
by 

(9 T/9 P)^ = l/(p S) (4.1) 

Afterwards, thermal energy is picked up isobarically. From the point of view of classical heat 
exchanger operation, isobaric assumptions apparently do not present very challenging tasks. 
However, the T-profiles of section 3 provide hints about design constraints, different case. 
Further, there is nearly no quantification at all for the real depressurization process at d|i = 0. 

The effectiveness results appear to throw some light on the refrigeration process. Ini- 
tial motivation for the vortex refrigerator introduced by Staas and Severijns (compare also 
App. A), appears to have been related to attainment of low temperatures. Obviously, the 
theroretical calculations based on ideal liquid do omit the lattice contributions of the real 


16 



"superleak" plug which provides the cooling. The refrigeration capability is elucidated by a 
brief look at the opposite direction i.e. a power cycle. Very little power is gained as the area 
enclosed by the cycle in a pressure-volume diagram is extremely small. This feature relates 
to the small volume changes in liquids. Thus, there are very few attractive points for power 
cycle operation. However the opposite appears to be true for refrigerator operation. Very lit- 
tle work is needed for heat pumping in view of the small volume changes. Moreover, the use 
of parasitic heat is a very attractive feature. It is near at hand to omit system preparation for 
a first look. Thus, the usual environmental temperature of about 300 K has been disregarded 
completely, (as in similar analysis of the fountain effect pump thermodynamics). The 
definition used for the COP (= coefficient of performance) is the ratio of heat pumped to heat 
input. This COP ideally is 

COP = T S /Q r (4.2) 

The simplest case is a refrigeration load 

Qr = T 4 S 4 + IQuI (4.3) 

Assumed is the sequence 1-2 for d|i = 0 and so on. The system’s thermodynamic path is 
shown schematically in the T-S diagram. Fig. 4.2. The coefficient of performance (COP) is 
depicted in Fig. 4.3. It is seen that reasonable values of COP are reached at relatively high 
bath temperature T b for a specified pressure difference. Fig. 4.4 shows the cold end tempera- 
ture versus T 2 . 

The effect of a counterflow heat exchanger (Fig. 4.5) is known from classical systems. 
The COP usually is lowered, if all other parameters are kept constant, and the heat exchanger 
covers a specified temperature difference prior to the refrigeration stage. The final tempera- 
ture is decreased for such a system, and the compensation of a drop in COP is possible in pri- 
ciple using a lower AP in the FEP. Figure 4.5 presents designations for the cycle with 
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Fig. 4.3. COEFFICIENT OF PERFORMANCE (COP). /D OF 

VORTEX REFRIGERATOR VERSUS INITIAL TEMPERATURE T, 
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Fig. 4.4. 


Ideal cold end temperature T 4 achieved in the p.-P-p.-P VR 
cycle, versus initial temperature Tj ("bath temperature"), for 
various pressure increases in the fountain effect pump (ideal 



precooling by counterflow heat exchange. Fig. 4.6 depicts the system. Fig. 4.7. shows the 
COP results and cold end temperatures for a parametric study with the precooling temperature 
ratio as parameter. A relatively small effect of precooling appears to provide considerable 
modifications of the ideal system performance. The parametric studies have been conducted 
for a pressure difference of 100 milli-bar. In the limit of a temperature ratio of unity, the pre- 
vious simple cycle |i-P-|i-P is recovered. 

Real component selection has to consider two items in each plug for ideally constant 
chemical potential. One case is the porous media selection for the specifications. The other 
case is the heat absorption/heat rejection unit optimization. In addition there are connecting 
plumbing devices and controls. The simple wide tube connection appears to be a low flow 
resistance solution. As shown by the FEP duct data however, there is a good chance of heat 
flow by counterflow action. Fortunately, there is a very small slip of normal fluid and 
superfluid at high mass flow Reynolds numbers. This implies a possible need for balancing of 
slip versus pressure drop in a real refregerator. Another solution is the use of "pressure 
relieved mass transport" using porous media coatings (wick-like structures) for wall coverage. 
There appears to be very little information available in this area. 

In summary, the vortex refrigerator has an attractive potential for He II - cooled sys- 
tems and subsystems, however considerable R & D work is required for general acceptability. 
This includes the interfacing between He I portions and He II subsystems during system 
preparation. Nevertheless, there are very advantageous features indicated by the present cal- 
culations which pint out good performance when T is not too low, e.g. near 1.8 K to 2 K. 
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FIGURE 4.5. Schematic diagram of cycle with counterflow 
HEAT EXCHANGER (HE) a. T-S diagram, schematically, 

log-log coordinates; 

b. log P versus T diagram, 

schematically; point 1 is on the 
saturated vapor-liquid equilibrium 
curve. 






Y 


CB 


Fig, 4.6 . VORTEX REFRIGERATOR WITH COUNTERFLOW HEAT 
EXCHANGER (CFHE) , schematically; 

AC AFTERCOOLER 
CB COLD BOX 

ESL "EXPANSION SUPERLEAK" = MECHANO-CALORIC UNIT 
FEP FOUNTAIN EFFECT PUMP WITH HEATER ACTIVATION 
V VALVE 
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Fig. 4.7. COP FOR CYCLE WITH IDEAL HEAT EXCHANGE ; 
INSET : COLD END TEMPERATURE VS T ] . 
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5. POROUS MEDIA TRANSPORT CHARACTERISTICS 


The experiments have included basic studies of vortex shedding processes during heat 
flow in porous media. In conjunction with the power law fits for FEP pump plug data, the 

determination of the rate of mass/superflow versus applied temperature difference has been a 
major goal of the present work. However, from literature the most common case of Zero Net 
Mass Row (ZNMF) been found to be pursued only in a limited range. In particular the large 
AT-range (non-linear range) has received but little attention. Therefore, this part of the stu- 
dies has aimed at an extension of heat transport knowledge, available previously from VLPS 
work. 


The experiments have used plugs whose axis has been oriented in vertical direction. 

Heat flow has been directed in the "g" direction downward. Heating has been initiated using 

step inputs in heating power. Heater windings have been mounted in an insulated chamber. 

The step input in power in general produces the following events: 

1. At low power, a nearly exponential first order change with a temperature difference AT 
proportional to [1 - exp (-t/x)] is observed; x = characteristic time; t = time. 

2. At a certain power, linearity is destroyed. The final value observed at long time is 
considered the result of quasi-steady operation. The q versus final AT curves are non- 
linear. 

3. Above a moderate power, there is no chance to reach a well-defined quasi-steady heat- 
ing condition at long times. Instead, at an intermediate time, there is a drop in AT as 
a result of vapor formation. 

4. Near the onset of case 3, there is a finite chance for large amplitude, slow, non-linear 
oscillations. 
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Figure 5.1 is a schematic drawing of the setup. Figure 5.2. presents a thermogram 
example. Figure 5.3 is a similar example with oscillations. Fig. 5.4. displays the results, at 
long times of quasi-steady linear and non-linear conditions. 

Figure 5.4 shows a linear region given by the permeability of the plug. The analog of 
Darcy’s law is used. 

qo = p S T v no (5.1) 

with a superficial normal fluid velocity of 

V* - I v PtVH. (5.2) 

The normal fluid permeability is k„. The grad P T - term in Equation (5.2) is the London pres- 
sure gradient of the thermomechanics of He II. 

Figure 5.5. shows permeabilities obtained as a function of temperature for plug PK -K-10-S- 
02- 6.4 x 0.75 

The interpretation of the linear regime is straightforward noting that the small AT of that 
range does not present difficulties of property definitions. 

In the non-linear regime things are less straight forward. A procedure has been used 
initially in parallel to phase separator (VLPS) work. The simple modified plug equation 
incorporaties a size-dependent (dimensionless) plug constant Kznmf* The equation is semi- 
empirical because of the implications of power law exponents discussed above in the preced- 
ing section. The simple analog of VLPS is given as superficial heat flux density q. 

q = k znmf ps S T f(p s /p n ) I V P T I (q n /p)] 1/3 (5.3) 

Results for the plug PK-K-10-S-02-6.4 x 0.75 are shown in Figure 5.6. 

Details of related numerical work are given in the Appendix Sections D and E. 
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Schematic diagram of porous plug 

Th- Thermometer, H- Heater, PP- porous plug 
P- phenolic, FT- Feed Through 

I - Insulation 


Fig. 5.1. Schematic diagram of experimental system for 
zero net mass flow (ZNMF) ; plug of PK series. 
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THERMOMETER VOLTAGE 

ARB. UNITS 



TIME, ARB. UNITS 


Fig. 5.2 THERMOGRAM EXAMPLE : PARAMETER : HEATER CURRENT 
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Eqs. MOTION 

- 2-FLUID 

MODEL , 

truncated, 


WILKS - 

BETTS , 

op- cit. 

^V s /bt = 

S grad T 

- <i/p) 

grad P 

c)v /^t = 
n 

- ( P '*'?*> 

S grad 

T - (1/p) grad P 


Nonlinear region 


In Linear region 

Slope = S 2 TX n s 46 - 42 mv, / cm 

~X* Z 

= 2.06 10" 10 cm 2 


FIGURE 5.4 Heater power q vs AT at bath temperature 
near 1.9 K. 

Inset : TWO-FLUID EQUATIONS . 




X. 

5x10 


1 xl 0 


5x10 


0 



Fig. 5.5. DARCY PERtCABILITIES VERSUS T . 
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Fig. 5.6. 


Rate constant 



as a function of bath temperature. . 


32 



The change to a more meaningful power law exponent involves a characteristic length. 
For instance, the case of an exponent of 1/4 is quite significant. There is a characteristic 
length involved in the interaction of pore walls with fluid motion. Appendix E presents 

details of this approach. In the course of the work, the data obtained indicate that the con- 
stants replacing in Equation 5.3 retain the size dependence qualitatively. 

In phase separator (VLPS) plugs the size dependence is at variance with ZNMF in so 
far as data appear to be numerically above the VLPS plug equation. However the data trends 
of the constants of Eq. 5.3 are retained. A comparison is available using Figures 5.7 and 5.8. 

In summary, it is seen that early indications of transport phenomena in plugs are 
confirmed. Jeffrey M. Lee had already indicated in his M.S. thesis (UCLA 1983) that certain 
plugs pointed in this direction. The VLPS results of Yuan’s Ph.D. thesis (UCLA 1985) put 
this suspicion on a firm footing in the VLPS area. The data obtained so far for ZNMF appear 
to confirm the trend of the VLPS data also for ZNMF. Thus, the peculiarities of FEP units 
are in considerable contrast to the conditions of VLPS and ZNMF modes. 
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Fig. 5.7 POROUS MEDIA RATE CONSTANTS OF PHASE SEPARATOR PLUGS VERSUS PERMEABILITY 




6. SUMMARY INCLUDING COUNTERFLOW COOLING OF COMPONENTS 

Studies of fountain effect-induced phenomena and system options have been continued 
with the goal of improved fluid management starting near 2 K. Instead of looking at fluid 
management only from the point of view of routing "hot" to "cold" production lines, the alter- 
nate scheme of managing from "bottom" to "top" is suggested by the He II work involving 
local He I domains. For instance, solid hydrogen pellets may be carried hydraulically instead 
of pneumatically on the basis of the Low-T cryogenic fluid. 

The porous plug results are an important data set for various thermomechanical modes. 
A temperature difference is not possible unless the porous medium restricts the cross section 
to suppress the flow of normal fluid. The permeability range covered in the present runs has 
been reduced from the previous phase separator range to the order 10 -10 cm 2 . The phase 
separator modes of previous work have been compared to the data sets for zero net mass flow 
(ZNMF). The ZNMF data do indeed follow phase separator plug trends. The absolute 
ZNMF transport rates obtained however are quantitatively higher than the phase separator 
data. 


The extension of the work to He II vortex refrigeration has been based on the preced- 
ing calculations of FEP units. The addition of the mechano-caloric plugs presents a novel 
component set to be quantified in transport details for the purposes of refrigeration and heat 
pumping. Previous vortex refrigerator work appears to have been motivated primarily by 
attainment of the lowest possible temperatures. The performance figures appear to be a spe- 
cial class among various cryocooler systems. The COP values are attractive for operation 
near 2 K i.e. in the upper T range of He II, and not in the lower one. 

Finally, there is an interesting implication for "lifetime extension". The usual space 
vessel of the IRAS type has the VLPS location set rather far away from heat input locations. 
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This approach appears to be advisable for the understanding of the mode of operation prior to 
more sophistication. With the success of IRAS however, the question for improvements 
includes the use of the VLPS effect and TM effect in order to intercept certain heat inputs. If 
less heat arrives at the bath, the vaporization rate is reduced and the He II is kept in the 
vessel for a longer period of time. Another point is the incorporation of the vortex refrigera- 
tor into this type of concept for improved fluid management of the future. 
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APPENDIX B: Vortex Shedding Regime in Porous Media: 

Modified Mutual Friction Transport Regime of Counterflow 

It is noted that porous media transport is not the usual mutual friction transport pro- 
posed by Gorter-Mellink. However, for an assessment of the special porous plug conditions, 
the frame of reference of Gorter-Mellink provides insight and empirical equations for VLPS 
design. First, vapor liquid phase separation (VLPS) is discussed. Second, zero net mass flow 
(ZNMF) is treated. 


Consider VLPS modes in the form proposed by S. W. K. Yuan. The rate constant K* 
* Kvlps is contained in the equation for the superficial mass flux density (= mass flow rate 
per total plug cross section divided by the density). 

To = Kv LPS [^+ST]- 1 p s ST[(p s /p)S(dT/dz](Ti n /p)] 1/3 

(B.l 

In general, Eq. (B.l) has to be integrated in order to take the T-dependence of properties into 
account. Details have been given in the Ph.D. thesis of Dr. Sidney Yuan. 


The related equation for the heat flux density associated with the mass flux density is 
given by q 0 = Xj 0 , or 


(q 0 )vLPS = K^ LPS X[X+ST]- 1 p s ST[SIVTI(p s /p n )0l n /p)] 1/3 

For finite temperature differences, Eq. (B.2) is integrated for steady transport: Jq^ dr. = 
KvLPsJf(T)dT, or 


(3o)vlps - K VLPs(L) 1/3 j 


J>. 3 a+ST) 3 (p s ST) 3 S(p s /p n )(Ti n /p)dT 


1/3 


(B.3) 


For the linear regime (Appendix T ), the relationship between VLPS and ZNMF is: 


(Qo)vlps - (Q 0 )znmf ,/ ( 1 +ST/A.) 


(B.4) 


PRECEDING PAGE BLANK NOT FILMED 

3 - '-fcT 
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Using Eq. ( " .4) for ZNMF results in the following: 

©ZNMF = KiMFP«ST[(VP)(P/Pn)(S AT/Az )] m AT « T 

(B.5) 

(Az = L). 


Again the case of large AT necessitates integration due to the strong temperature 
dependence of the properties p s (T) and S(T). By grouping all temperature dependent proper- 
ties into one term F(T), Eq. (B.5) becomes: 

CUZNMF = Ki, MF [F(T) vn lfl 
The property function is 


F(T) = p s Wry(P„P> 


Integration leads to 

zs=L T+AT 

J (Q 0 )znmf^ z = k znmf J F(T)dT 

z=0 T 

The result is 


(B.7) 


(B.8) 


GWznmf - (K*) ZNMf L- w [jF(T)dT] 1 ' 3 
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APPENDIX C: EXPERIMENTAL DATA FOR PLUG PK-K-10-S-02-6.4 x 0.75 


(raw data) 


The tables contain the temperature differences across the plug 
and the superficial heat flux density (heat flow rate per 
total plug cross section). 


dT,K 

0.011 

0.0198 
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0.096 
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q,mW/cm 2 
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AT,K 

q.mW/cm 

0.00125 

34.30000 

0.00120 

32.04000 

0.00150 

56.97000 

0.00156 

81.90000 

0.00200 
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4 T,K 


0.0005 

78 . 1700 

0.0009 

94.1700 

0.0010 

104.9000 

0.0016 

125.7400 

0.0018 

139.4000 

0.0026 

149.7700 

0.0036 

151.1000 

0.0046 

157.7800 

0.0086 

167.4100 

0.0008 

173.0400 

0.0111 

193.4600 

0.0141 

215.0000 

0.0161 

227.8600 

0.0236 

244.4300 

0.0266 

247.8100 


T = 2.00 K 





Q.CO0 


0.012 


D.01B 


n.02ti 


4T,K 
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AT.K 

q,mW/cm^ 

0.0117 

116.2600 

0.0170 

149.7700 

0.0200 

181.6500 

0.0280 

219.8000 

0.0410 

261.5800 

0.0480 

279.3100 

0.0600 

297.6200 

0.1050 

356.0400 

0.1070 

376.6700 


q, mW/cm 


2 
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4T,K 

qtmW/cm 2 
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45.9200 
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65.3900 

0.0055 
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0.0065 
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0.0110 

122.1400 
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0.0165 

144.5300 

0.0255 

164.6300 

0.0310 

178.6700 

0.0450 
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231.1300 
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0.1300 

314.6000 
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AT,K 

q, mW/cm 

0.0010 

47.6100 

0.0020 

58.6100 

0.0020 

96.2700 

0.0010 

68.0300 

0.0015 

17.2300 

0.0022 
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106.2700 

0.0025 

116.2500 

0.0030 

100.5500 

0.0055 

164.6300 

0.0120 

191.9700 

0.0150 

216.6100 

0.0165 

216.6100 

0.0295 

254.6500 

0.0245 

242.7500 

0.0430 

277.5100 

0.0460 

282.9200 

0.0570 

308.8900 

0.0575 

308.8900 

0.0695 

320.3600 


q, mW/cm 


BATH 1.91 K 
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o 
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0.0245 
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244.4300 
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0.1575 

372.5000 


q, mW/cm 


BATH : 1.71 K 


400 


3S0 


300 


250 


200 


ISO 


100 


50 H 



AT, K 


51 




APPENDIX D: PERMEABILITY DETERMINATION 
BASIC EQUATIONS 

The permeability for ZNMF modes in the laminar region (i.e. low counterflow) is expressed 
within the generalized Darcy low as: 

V no = I V P T l/Hn (D.l) 

V no = Initial normal fluid velocity 
P T = Thermomechanical pressure 
Ti n = Viscosity of normal fluid 
Kn = Permeability 

London’s equation predicts that the theimomechanical pressure 
gradient can be produced externally by a temperature gradient: 

VP = pS V T (D.2) 

p = Fluid density 

S = Entropy 

VT = Externally applied temperature gradient 

Finally, for ZNMF modes, the heat flux density (q„) in 
the axial direction of flow through porous plugs is given by 

flo = P STv no 

so that through equations fD.l) and ( D.2), the governing 
equation for permeability calculations in ZNMF is presented 
as: 

q 0 = pST tc„ pS (AT/L)/n n 

where often AT «T, and there is a property function 
f(T) = ? 2 S 2 TM n 


(D.3) 

(0.4) 

(D.5) 
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APPENDIX E 


THERMOMETER CALIBRATION AT VERY SMALL 
TEMPERATURE DIFFERENCES 


The carbon resistor thermometers used are Allen - Bradley, 39 Ohm, resistors. The 
calibration procedure involves the measuring of the thermometer’s resistance and the 
corresponding vapor pressure of the system. The properties of the thermometer are such that 
the resistance is approximately logarithmically proportional to the vapor pressure. This can be 
expressed as: 

In (R) = a In (Pv) + b ( E . 1 ) 

R = Thermometer resistance 

Pv = Vapor pressure 

a,b = Linear least square parameters 

Therefore, by using a temperature-vapor pressure curve, such as the T-58 scale, the 
temperature versus resistance relationship is derived for small temperature difference. A prin- 
tout of a sample computer program and its output shows the utilization of a T-58 cubic spline 
fit (fit presented by Donnelly, R.J., et al in J. Low Temp. Phys., 1981) for iterative resistance- 
temperature calculations 

For very small T-differences, the calibration equation ( .1) Permits a convenient deter- 
mination of AT from the resistance difference AR. The slope of the vapor pressure curve is 
known from the T-58 scale: 

T = P v “ 1 dP v /dT < E - 2) 

Inserting d In R/d In Pv = a = (dR/R)/(dP v /Pv) one obtains 

AT = AR dT/dR = AR (aRT) -1 

AT « T (E. 3) 

R = R 0 


R = R 0 at the bath temperature. 
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APP. D : <contin. > , 

SUBJECT : CUBIC SPLINE FITS TO 

THERMOPHYSICAL FUNCTIONS 

Example s ABOVE THE LAMBDA TRANSITION 

Lit. i C. F. Barenghi, P. G. J. Lucas, 
R. J. Donnelly, J. Low Temp. 
Phys. 44, 491, 1981 . 

Selected smooth spline fit resultsi 

Spec. heat at vapor pressure of saturated 


liquid , 4 internal knots, 
scale; C sat in 

basis i T— 58 
J/ (mol K) 

* 

e* 

1 

1- 

T, K 

Csat 

10E-4 

2. 1721 

45. 777 

10E-3 

2. 173 

33. 537 

10E-2 

2. 182 

21.222 


The corresponding Prandtl numbers (Pr) 
are 0.2891 , 0.5487, 0.7685 . 


T-58 SPLINE FIT 


C.F. BARENGHI, R.J. DONNELLY, 

R.N. HILLS, J. LOW TEMP. PHYS. 51,319,1983. 
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APPENDIX F. ABSTRACTS OF RECENT CONFERENCES 


Am. Phys. Soc., Div. of Fluid Mechanics, 
40th Annual Meeting, Eugene Oregon, 22-24. 
Nov. 1987. 

Non-Linear Counterflow in Porous Media* 


HI 6 

Nonlinear Counterflow in Porous Media. P. K. KHANDHAR and 
T. H. K. FREDERKINO, University of California, Los Angeles, CA 
90024 - Sintered stainless steel plugs have been subjected to coumcrflow 
using liquid He II at zero net mass flow (ZNMF). The normal fluid flow, 
emitted from a heater, is compensated for by an inflow of superfluid into 
the heating chamber. The results give a permeability on the order of 
10" 10 cm 2 in the linear region based on a Darcy law analog (nominal 
filtration rating near 1 pm of classical fluid flow through the sintered 
plugs). The non-linear region established at increased heater power is 
characterized by a ZNMF rate constant K* (Gorter-Mellink mutual fric- 
tion constant). In contrast to the wide duct value, the present K* data are 
pore size dependent K* is found to be a monoton icaJly increasing func- 
tion of the permeability in qualitative agreement with porous plug data 
for vapor-liquid phase separation (VLPS). VLPS results for K* have 
been found to be proportional to the square root of permeability 1 . 

'S. W. K. Yuan. Ph.D. Thesis. UCLA, 1985. 


*) Supported in part by NASA AMES Research Center. 
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App. F (ABSTRACTS, con tin.) 


PERFORMANCE TEST OF A LABORATORY PUMP FOR LIQUID TRANSFER 

BASED ON THE FOUNTAIN EFFECT 

W.E.W. Chen, T.H.K. Frederidng and W.A. Hepler 

University of California, Los Angeles 
Los Angeles, Callifomia 


ABSTRACT 

A laboratory scale pump has been tested in detail in order to determine the flow 
characteristics of a heater-activated all-metal thermomechanical pump (fountain effect 
pump FEP). The emphasis is on the functional dependence of the fountain pressure 
difference versus mass throughput A modified Vote et al. power law approximation is 
employed for a simplified description of the flow rate as a function of the driving 
force. Flow rates of up to 10 Litres/(hr cm 2 ) have been obtained despite a large nomi- 
nal pore size of the porous plug of 2 pm (filtration rating) used for the FEP. 


1987 Cryog. Eng. Conf. St Charles, IL; June 1987, paper BC-6 
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APPENDIX G 


REFRIGERATION AND HEAT PUMP SYSTEMS BASED ON He II VORTEX CONTROL* 
T.H.K. Frederking, H.H.D. Tran and R.M.Carandang 
University of California, Los Angeles 90024 CA 


The utilization of thermomechanical forces permits 
pressure increases in pumps without moving parts. This topic 
has become an interesting option for liquid transfer at micro- 
gravity. In addition, the thermomechanical pump (fountain 
effect pump) may be used to energize refrigerators, substitu- 
ting for the usual mechanical pump. In the present thermody- 
namic studies, cyclic variations of the fluid equilibrium 
state are Investigated Including performance parameters, in 
contrast to attainment of low temperatures, the emphasis Is 
on heat pumping tasks at temperatures between 1.5 K and the 
lambda temperature. 


Key words: He II; thermodynamic cycles ; fountain effect; 
refrigerator/heat pump systems. 


I. Introduction 

Increased use of liquid He II as a coolant near 1.8 K has included application 
of the fountain effect (thermomechanical effect). Examples are vapor- liquid phase 
separators and fountain effect pumps (FEP). Another type of application addressed 
In the present studies Is the use of FEP - induced pressurization for refrigeration 
purposes and heat pumping. An original version of such a system, based on the con- 
trol of the dynamics of quantized vortices, has been proposed and tested by Staas 
and Severijns [l] L This Staas-Sever I jns cooler has aimed at the attainment of low 
temperatures down to 0.7 K . In contrast, the present work emphasizes heat pumping 
In the range of the temperature (T) encountered in space cryogenics and in supercon- 
ducting magnets for high magnetic fields based on NbTi/Cu. The system applicability 


Numbers in brackets refer to the I iterature references listed at the end of this 
paper. 

*) Presented at 4th Intern. Cryocooler Conf. Easton, Md. 1986. 
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Therefore, heat is to be removed in order to keep temperatures, at a specified loca- 
tion, time - independent. In the aftercooler (AC), the temperature is brought back 
to the bath temperature. In the "expander", application of the pressure difference 
generated in the FEP, forces superfluid through the plug. This results in cooling 
by the mechano-ca loric effect. The refrigeration load (Op) is supplied in the cold 
box (CB). Finally, the capillary (CAR connects the cold box to the bath. 

Figure 2 indicates a scheme with a valve (V) replacing the capillary system. 
In principle, other means aFe available, e.g. a valve at the inlet. This option is 
not favored for a saturated liquid bath, as it may lead to pressure reduction and 
cavitation. This danger is avoided with pressurized liquid He 1 1 . It Is noted that 
in Reference [I] a set of capillaries has been employed in order to test the flow 
conditions in the entropy carrying tubes downstream of the cold box. 

4. Fountain Effect Pump Performance 

The fountain effect unit is the driver for the vortex refrigerator. Therefore, 
its performance exerts an important influence on the coefficient of performance 
(COP) of the refrigeration process. Recent work on FEP use for liquid transfer in 
space [3] has prompted a critical inspection of various performance measures . 

Two types of effectiveness or efficiency figures have become known. One relates to 
the liquid yield, ultimately available at the receiver, during transfer from one 
vessel to another one. The other measure is the energetic effectiveness e e . Because 
of generally small values obtained for e , there has been a sufficient incentive 
to Improve the energetics of the FEP , e . g. by heat pump - assisted FEP units [4], 
Ideal changes of state at dp=0 are referred to in the subsequent discussion of the 
FEP. 

The energetic effectiveness is defined^as the ratio of the flow power at the 
pump outlet , Wp , to the heat Input rate , Q , energizing the pump. Neglecting the 



Figure I. Vortex refrigerator, schematically; AC aftercooler, CAP cap i I lary, CB cold 
box , ESL "expansion super leak", FEP fountain effect pump, H heater ; 
Refrigeration load Qp is absorbed in the cold box. 
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( 2 ) 


grad ( \i + v s ^/2) = 0 

where the kinetic term is small outside the narrow passages formed by the porous 
medium. The momentum equation of the model , = p v = p s v s + p n v n , for normal 

fluid locked in the plug, reduces to the simple condition 

P* =p s t 5 (3) 

Thus, the superfluid critical speed v sc is related to the critical velocity of the 
observed mass flow by v c = v sc p s /p 

In the terrestrial gravity field, the fountain effect is quantified readily by 
Inclusion of the gravitational contribution to the generalized chemical potential 
function ( p + gaz + v s ^/2 ) ; (z vertical position coordinate). In a thermo - 

static device. Equations (I) and (2) have in common the limit of the thermomechani- 
cal gradient, at zero flow, of 

grad P = grad Py = p S grad T (4) 

(London pressure gradient). When a jet geometry allows a speed v , downstream of 
the superleak, the limit of ideal flow without losses leads to a fountain height 
of Az = v* /(2 g) . 

There are two types of fluid processing components in general: First, the po- 
rous plugs needed for the thermomechanics and second, ducts and conventional de - 
vices, such as heat exchangers. In the latter category, significant vortex shedding 
may take place causing dissipation of flow power. At sufficiently high flow rates, 
both fluids may move with nearly the same speed. In this case, interaction of the 
normal fluid with the vortices causes "near - classical " flow. For Instance, in a 
duct with a diameter beyond the order of I0“2 cm, roughness - dependent friction 
may be observed readily, as in Newtonian turbulent fluid flow . Therefore, classi- 
cal equations have been used in the analysis of this aspect of the Staas-Severi jns 
refrigerator [l] . Alternate use of vortex pinning components and of vortex shed - 
ding components implies controlled processing of vortex systems giving rise to the 
designation of "He II vortex refrigerator" , ( or of " superfluid vortex fridge"). 
The various dissipation aspects in wide and narrow passages have been reviewed by 
Hammel [2] giving details of basic thermophysics of the He II modes of motion. 


3. Simplified Vortex Refrigerator Operation: Concepts 

In the discussion of this section, simplified refrigeration concepts are in- 
voked. Up to this time, little has become known about performance figures. In line 
with applied thermodynamics conventions, the coefficient of performance is defined 
as the ratio of the refrigerat ion load to the net work. The thermal efficiency of 
thermal to mechanical energy conversion Is defined as the ratio of the net work to 
the heat input. Macro-thermodynamics Is used. For Instance, at a superficial glance 
one may refer to the FEP device as "compressor". The reverse case of mechano-ca lo- 
ric cooling in a plug may be referred to as "expansion device", (or "ideal superleak 
expander"). Caution is needed though to avoid a literal interpretation In terms of 
classical ideas which may lead to non-physical consequences . 

A schematic diagram of the original system of the Staas-Severi jns vortex re- 
frigerator is shown In Figure I. There are "compressor" and "expander " peculiari- 
ties elucidated by the two-fluid model for du = 0 . "Compression" is actuated by 
the heater (H) which supplies thermal energy. As superfluid passes through the 
device at dp = 0, normal fluid left behind tends to accumulate entropy upstream. 
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appears to be not restricted to these examples and general features are presented: 
First, the vortex dynamics is discussed. Subsequently, simplified refrigerator oper- 
ation is considered. The next subject is the FEP pressurization device. Heat pumping 
performance is considered for one of the simplest cycle options, and an outlook is 
given. 


2. Vortex Dynamics and Two - Fluid Aspects 

The quantum liquid He II has a high heat transport rate in conventional refrig- 
erator ducts, unless the flow cross section Is severely restricted. Therefore T - 
gradients are small in usual flow passages. In order to generate sizable temperature 
differences, needed for the application of thermomechan ica I forces, fine porous 
media, called superleaks or superfilters, are used. The limiting case is the "ideal 
superleak" with zero heat flow. The porous medium permits res i stance I ess flow (*"su- 
perflow") up to Its critical velocity (v__). As the pore size Is reduced, at a 
specified temperature, the v - value attainable is Increased. At v sc , liquid 
circulation and flow resistance is Initiated. Because of quantization of the liquid 
circulation, a finite energy is needed for vortex motion onset with velocity com - 
ponents perpendicular to the main flow direction, I. e. vortex shedding . The energy 
needed is large for very fine pores. For FEP purposes, critical velocities of the 
order of 10 cm/s have been considered convenient In laboratory work. Because of the 
energetic restr ict Ions, the fine plugs may be regarded as "vortex pinning devices" 
below v sc . In the ideal thermodynamic treatment of changes of state, the ideal 
superleak constitutes a perfect vortex pinning component. 

Once a vortex Is activated above v sc , a Magnus force provides vortex shedding 
e.g. from walls toward the fluid Interior. This phenomenon leads to a finite differ- 
ence In chemical potential (p) *. In the Ideal super leak dp * 0 is observed up to 
the critical velocity. Small vortex shedding rates of real plugs are expected to 
cause only small departures from the Ideal thermodynamic change of state. Therefore, 
ideal processes involving changes of state at constant chemical potential are use- 
ful for performance assessment. It appears that minor vortex shedding rates cause 
only small departures from the ideal performance values of FEP units. 

We have a flow system described In usual hydrodynamic terms supplemented by 
the superflow thermodynamics. The two - fluid model Is useful . The existence of 
two interpenetrating fluid media is postulated: There exists a viscous normal fluid 
component and an inviscid superfluid component ( a+"concen+ra+ions " ( p n /p ) ar| d 
(ps/'p) adding up to unity). The flow field Is subdivided accordingly: There is a 
normal fluid flow pattern and a superfluid flow field. Because of the fine pores, 
the normal fluid motion is considered negligible, in particular in the limit of 
perfect immobilization of superleaks. Thus, in the present context , attention is on 
•the superfluid described by 

Dv s /Dt = -Vp *7P/p - S7T - (\/2)(p n /p) Iwf N) 

with w = v n - v s = relative velocity ; S entropy per unit mass; p liquid density). 
For steady, one-dimensional flow with negligible |w£ term, we have a simple truncat- 
ed form of the preceding equation 
_ 

In part of the low temperature I iterature, the Gibbs free energy is customarily 
denoted as p . In the present discussion , this notation is adopted with p chem- 
ical potential per unit mass. 
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'Figure 2. Vortex refrigerator, schematically; V valve, (other components as 
Figure I). 


in 


mass loss and using Equation (4), one obtains 

e e = W F / Q = SLT /ST (5) 

S Is the mean entropy for the range from the bath temperature T to the temperature 
(T + AT) at the pump exit. For very small temperature differences, S becomes equal 
to the bath value S(T), and the effectiveness simplifies to 

e g = AT/T ; AT <C T (6) 

A comparison with the thermal efficiency of a Carnot cycle shows that the FEP ef - 
fectiveness, in this limit, is equal to that efficiency for a low temperature en - 
vlronment. More specifically, the comparison Carnot power cycle receives heat at 
the temperature (T +AT) and rejects heat at T. 

Figure 3 presents a comparison of the FEP effectiveness with the thermal effi- 
ciency (qj) of the related Carnot power cycle. According to Carnot's theorem, the 
«W - values constitute the best performance results available Ideally. It is seen 
that E e is always below f|j . The effectiveness Is about 50 % lower than the Carnot 
limit. There is a common asymptote at low AT. The parameter adopted in Figure 3 
Is the bath temperature upstream of the FEP unit. At sufficiently high AT, the FEP 
performance tends toward an asymptote, near 15 % of the ideal heat pump-assisted 
FEP unit [4] . 

Figure 4 shows the heater power requirements. Figure 4a is a plot of the heat 
Input rate per mass flow rate (Q/m), starting at a temperature of I .6 K versus the 
temperature difference across the ideal pump. The upper curve Is the Ideal require- 
ment of the simple FEP . The real pump will need a somewhat larger power because of 
departures from the ideal change of state characterized by dp = 0 . The lower 
curve depicts the power requ irements, W/m , per unit mass, of the heat pump-assisted 
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Fig. 3. Energetic effectiveness^ and thermal efficiency ru. of Ideal systems for 
various bath temperatures; (HPA-ISL heat pump-assIsTed Ideal superleak). 
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Fig. 4. Heater Input (Q/m) and power input (W/m) functions for ideal FEP systems; 

a. power required versus £T , bath temperature 1.6 K; b. parametric plot of 
the power required vs. T for various bath temperatures. 







FEP unit. This power is low at low and rises monoton i ca I I y with an increase in 
the temperature difference. 

Figure 4b summarizes the power input requirements in one graph making use of 
a parametric representation with the bath temperature (upstream temperature) as pa- 
rameter. It is seen that there is a unique Ideal heater power function (Q/m) given 
by the ideal superleak constraint : (Q/m) is the value of S T at the pump exit. For 
instance, for a mass flow rate of m = 40 g/s , or a volumetric flow rate of about 
1000 liters/hour, a heater power of the order of 40 Watts is required near 1.8 K. 

At large temperature differences, the heat pump assistance is no longer attractive 
as the power requirement approaches that of the simple FEP unit. 


5. Refr igeration - Heat Pump Cycle Example 

One important assumption in the preceding and in the present treatment is at - 
tainment of local thermodynamic equilibrium . Other system conditions of the present 
section are as follows: steady operation, steady flow of liquid and of the two 
fluids, simplified changes of state using ideal conditions, such as isobaric changes 
in heat exchangers, and ideal superleak processes at constant chemical potential. 

A simplifying tool is the introduction of the Ideal heat pump-assisted FEP (dp= 
0) and a related state change In the "expander " subsystem. Furthermore, for quick 
estimates the contribution from the integral / S dT may be neglected in the T-range 
near 1.8 K, as discussed elsewhere [4] . The changes of state are listed in Table 1 


TABLE 1 . STATE CHANGES OF THE SIMPLE CYCLE p-P-p-P 


juence 

Constra int 

' “ j 


1 - 2 

dp = 0 

2-3 

dP = 0 

3-4 

dp = 0 

4 - 1 

dP = 0 


Subsystem/Component 

Ideal FEP, heat pump-assisted 
isobaric aftercooler (T^ = Tj) 

Ideal superleak"expander" 

Absorption of the refrigeration load Q R 


The cycle is sketched in the T-S plane and in the P-T plane In the Insets of Figure 
5. After processing of fluid in the FEP, the aftercooler brings the temperature 
back to the bath temperature Tj, Subsequent "expansion" brings about a low temper- 
ature, and finally the change of state from T 4 to T| Is assumed to be available for 
the absorption of the refrigeration load. 

The temperature T 4 has been plotted for various pressure Increases in Figure 5. 
For instance, the pressure rise in the FEP unit, starting from 1.8 K, produces a 
T 4 -value near 1 . 6 3 K; (pressure difference achieved across the FEP of about 100 
mi lli-bar). Figure 5 includes the equivalent Carnot cycle temperature, for an "en- 
vironment" of 1.8 K, which has the same ref rigeration load as the p-P-p-P cycle 
under consideration. Both temperatures are displayed versus Q R . 

Coefficients of performance of the Carnot cycle (COP)cARNOT are shown in Fig. 6 
versus T . As the initial "environmental temperature" , e.g. I .8 K , is lowered by 
a small reduct ion AT T, the COP is the reciprocal qj - value of the power process. 
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Fig. 5. Lowest temperature of the Ideal p-P-p-P cycle and equivalent Carnot cycle 
temperature T c vs. refrigeration load Qr/iti ; "environmental" T * 1.8 K = 
initial bath temperature Tj j Inset : Cycle in the T-S diagram, schematically 
and cycle in the P-T diagram , schematically. 
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Fig. 6. Coefficients of performance of the Carnot cycle, (COP)c/\rnqt versus T for 
d i f ferenf'environmenta I temperatures* Te . 






6. Perspectives 

Simplified changes of state, with heat pump assistance of the FEP process, and 
a similar "expander assisted" pressure reduction, have been considered. The imple- 
mentation of near - ideal conditions by external means may require considerable com- 
ponent developments. Therefore other approaches toward improvement of the very 
simple cycle ought to be considered. One possibility Is a counterflow heat exchanger- 
Another one Is via"Carnotization" by means of multi-stage operation. 

Figure 7 depicts a schematic diagram of a vortex refrigerator with a counter- 
flow heat exchanger Increasing the T -range covered by the heat pump process. 

The counterflow avoids excessive entropy production brought about otherwise by a 
large temperature difference between the "high pressure" stream and the"low pressure* 
stream. 

The other option of Carnotizat Ion Is brought about by multi-stage "compression" 
and mu!ti-stage"expansion" . Figure 8 presents an example of a three-stage FEP - 
system with aftercoolers (AC). The inset of this figure is a sketch of the multi- 
stage cycle . As the number of stages is increased, the Carnot case is approached 
more and more. While the low density change makes multi-stage "expans ion"*less at- 
tractive, the possibility of multi-stage "compression" has received attention, e.g. 
Ref. [5] for the special case of He^ - He^ dilution refrigeration, and Reference [6] 
for the T-range near I .8 K In the space cryogenics work of recent years. 

In the final assessment, a major point is the FEP performance. The ideal ener- 
getic effectiveness figures of at most 15 % are not very attractive. A major advan- 
tage of the FEP - supported system is the general ly high reliability due to the ab- 
sence of moving components (aside from valve positioning). The most attractive 
issue appears to be the possibility of"waste heat"uti I ization (e.g. vessel heat 
leaks. Instrument dissipation rates). For high field magnets, an example of liquid 
circulation has been presented in Ref. (7]. For space cryogenics systems, the use of 
a cooler-activated FEP unit Is a viable option. Also heat might be"piped" from 
solar energy concentrators to the FEP without the use of photo-vo I ta Ics. Further 
recent results are in Ref. [8] . 

In summary, attractive features of the He II thermomechanics are available for 
specific refrigeration and heat pumping tasks. 
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Figure 7 . Vortex refrigerator with a counterflow heat exchanger (CFHE), schematical- 
ly; AC aftercooler , CB cold box ; ESL "expansion super leak" , FEP 
fountain effect pump with heater; V valve . 

*) from a "classical" point of view. 
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Figure 8. Example of multi-stage "compression" with three stages of FEP units de- 
signated as I, II, III ; AC after coolers In He II supply tank (schemati - 
cally ; Inset : Temperature - Entropy diagram for multi-stage "compress! 

on and "expansion" , schematically. 
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ADDITIONAL COMMENTS 


Numerous details of dissipative processes including narrow duct 
phenomena have been presented in the survey paper of E. F. Hammel 
at the Int. Conf. Low Temp. Phys. LT-10 , 1966 . This work has been 
based on the extensice studies of Keller and Hammel (also documented 
in the monograph of W.E. Keller on Helium-3 and Helium-4, Plenum Press 
1969. 

Concerning thermodynamic efficiency figures and similar measures, 

H.B. Callen*presents several examples : (simplified expressions on the 

right hand side are Carnot 

Thermodynamic engine efficiency values) 

Z e = d W /j-d Q h | = 1 - T c /T h 


Coefficient of refrigerator performance 

£ r = I- dQ c l/ I- d W | *T C / (T„ - T c ) 


Coefficient of heat pump performance 


t p = dQ / | - dW | 


T H / < T H ' V 


*) THERMODYNAMICS AND AN INTRODUCTION TO THERMOSTATISTICS, 1985, 
John Wiley, New York. 


67 



